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Abstract 



We obtain a general exact solution of the Einstein field equations for 
the anisotropic Bianchi type I universes filled with an exponential-potential 
scalar field and study their dynamics. It is shown, in agreement with previous 
studies, that for a wide range of initial conditions the late-time behaviour 
of the models is that of a power-law inflating FRW universe. This property, 
does not hold, in contrast, when some degree of inhomogeneity is introduced, 
as discussed in our following paper II. 



1 Introduction 



Two of the present authors have recently studied exact model universes filled 
with an exponential-potential scalar field and have found fO , a persistent 
non-inflationary behaviour in some exact Bianchi types III and VI, as well 
as in some cosmological models which may be thought as inhomogeneous 
generalizations of Bianchi type I cosmologies. These studies cast a doubt 
on the so-called "inflationary paradigm" which, while never properly formu- 
lated, vaguely states that the universe had undergone a period of inflationary 
expansion which is not only a must to solve a host of problems of standard 
cosmology, but is rather a "typical" dynamical behaviour common to a wide 
class of scalar field cosmologies. 

One of the main reasons the inflationary universe scenario is considered 
to be attractive || is that it has brought back to life Misner's hope [|J to 
explain the large-scale homogeneity, implied by the measurements of the cos- 
mic mirowave background radiation, without a need to impose very special 
conditions on the initial expansion of the universe. Strangely enough, how- 
ever, most of the work on inflation is done in the framework of isotropic and 
homogeneous FRW universes (Olive || and the references there in). 

This paper deals with exact solutions of the Einstein field equations. We 
show how one may obtain a general exact solution for Bianchi type I cos- 
mologies with an exponential-potential scalar field. Apart that one obviously 
needs a specific model for the potential to solve exactly the Einstein equa- 
tions, it is convenient to concentrate on this sort of the potential, often used 
in inflationary analysis, for the following main reasons: 

1. Exponential-potential scalar field introduces rather a small additional 
degree of nonlinearity into Einstein field equations, as it will be explained 
later, so that possibilities exist to solve these equations analytically in a 
variety of cases. In the case of cosmological models with Bianchi type I 
spatial symmetry one can obtain a general solution. 

2. It seems that the exponential potentials for scalar fields arise || in 
many theories such as Jordan-Brans-Dicke theory, the superstring theory, 
Salam-Sezgin theory and others. 

This work represents an introductory step towards our following paper 
II where the homogeneity will be broken in one direction to study the in- 
homogeneity effects on the late-time behaviour of scalar field cosmological 
models. 
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The main result of this paper is to confirm, on the basis of exact solutions 
we obtain, previous qualitative and numerical studies by other authors and 
the pattern of the late-time behaviour of the exponential-potential scalar field 
cosmologies. This agreement is only true, however, as long as the underlying 
geometry of the cosmological model is as simple as that of Bianchi type I. 
When the spatial symmetry group is more complicated, as for example in 
Bianchi III and Bianchi VI models as well as when the spatial homogeneity 
is broken the dynamics of the models is different. This, however, will be 
discussed in a separate paper II. 

In the following Section II we discuss the geometry and the matter content 
of the models. In Section III the Einstein equations are considered and solved. 
Some representative solutions are given in the section IV, while in the Section 
V the late-time behaviour of the generic solutions is studied qualitatively and 
numerically. In Section VI we conclude and summarise our results. 



2 The Bianchi I Exponential- Potential Scalar 
Field Cosmologies 

The usual synchronous form for a general Bianchi type I element is given by 

ds 2 = -dT 2 + a\(T) dx 2 + a\ (T) dy 2 + a 2 3 (T) dz 2 , ( 1 ) 

representing the anisotropic generalization of the spatially flat FRW uni- 
verse expanding differently in x,y, and z directions. While the form (1) is 
frequently used to study the dynamical behaviour of the cosmological model, 
we have found it convenient, in order to obtain exact analytic solutions, to 
cast the metric into a following, so-to-say, semi-conformal form 

ds 2 = e m (-dt 2 + dz 2 ) + G(t) (e pit) dx 2 + e~ p ^dy 2 ) . (2) 

This form of the metric is of particular interest when studying the inhomoge- 
neous generalizations of the line element (1), and is obtained from the Eq.(l) 
by the following identification 

dt = — , = a\ , G = a\a 2 , e p = — . (3) 
a 3 a 2 
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Since from now on we will be working with the line element (2) exclusively 
it is worthwhile to specify some cases of interest corresponding to particular 
values the metric functions / , G and p can take. 

When the function p vanishes globally the spacetime expands equally 
in axes x and y in which case the spacetime is called Locally Rotationally 
Symmetric (LRS). If in addition to p = one has G = e-^ then the metric 
is of the FRW form. In fact the only way the metric (2) may approach the 
FRW solution is when e-^ ~ G together with p ~ 0. 

The line element (2) admits important vacuum solutions 

K 2 - 1 

G = t, p = K\ogt, f= logt, (4) 

which are of Kasner type and always anisotropic apart from two "degenerate 
cases" k = ±1 which are disguised Minkowski line elements. 

The energy-momentum tensor for the scalar field driven by the potential 
V{4>) is given by 

T a(3 = 0, Q i/3 - g a p Q0, 7 0' 7 + V((f>)j , (5) 

and as long as the attention is concentrated on the homogeneous space-times 
the Eq.(5) may be rewritten in the perfect fluid form 



where 



together with 



T af 3 = ip + p)u a u p + pg al3 , (6) 



u a = - '= . (7) 



P = ~\<Pa^ + V(<f>) 

P = -l^-V^) (8) 

The perfect fluid interpretation of the scalar field, while not obligatory, is 
very useful to study the kinematical behaviour of the cosmological models. 
It is convenient to introduce the expansion 

e = u M; „<r (9) 
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and the deceleration parameter 



g = - 3 e 2 [e ]a u a + -e 2 J . (10) 

The sign of the deceleration parameter indicates as to whether the cos- 
mological model inflates. The positive sign corresponds to "standard" decel- 
erating models whereas the negative sign indicates inflation. 

For the models to isotropize at late times the shear measured with respect 
to expansion rate must vanish asymptotically . The models defined by the 
line element (2) are simple enough so that one may get a precise idea as 
to whether the spacetime isotropizes just from looking at the asymptotic 
behaviour of the line element. 



3 Solving the Einstein Equations 

The Einstein Equations for the metric given by the line element (2) and the 
matter specified by the stress-energy tensor (5) are given by 

* + g* + C ^ = 0, (11) 

§ = 2e^, (12) 
P+^P = 0, (13) 

!-H!) 2 -l /+ ? 2= -^ (14) 

where the Eq. (11) is the Klein-Gordon equation for the scalar field. 

Note that the Eq. (13) is somewhat decoupled from other equations in 
the sense that it is the same as in the vacuum case. The matter field acts 
indirectly on this equation through the "transitivity area" function G. Yet, 
one may immediately integrate the equation in terms of the function G: 

a 

P=q, (15) 
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where a is an arbitrary integration constant. 

Now we substitute the Eq. (15) into the Eq. (14) and use the fact that 
on differentiating the Eq. (12) and using it again one gets, for V ^ 0, 

f G G V 

f = G~G-V- (16) 
We then are left with the two following equations 

1 G 1 1 G dV 

* + G* = -2GV8t (17) 
G 1 G 2 + a 2 G G G V l2 

G + 2^-M + cr-^ (18) 

where the Eq. (17) has been obtained by substituting the Eq. (12) into the 
Eq. (11). 

We now specify the potential 

V{<f>) = Ae^ (19) 

Note that the Eqs. (17) and (18) are not valid for the case V = 0, in this 
case, however, the original system of equations is much simpler and its general 
solution is given by 

A = 0, G = t,p = a\ogt, <j> = piogt, f = ^ log*. (20) 

When a = and (3 = ±\/3/2 the above solution describes an isotropic 
spatially flat FRW universe. 

Returning now to the general case A ^ 0, the Klein-Gordon equation (17) 
takes the form 

v G ■ kG 

*+G*=-2G- < 21) 
In analogy with the Eq. (13) one readily obtains from the Eq. (21) 

m k G . 
*=G-2G' (22) 

where m is an arbitrary constant. 
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It is important to note that the exponential potential saves one much 
trouble with the non-linearity of the Klein-Gordon equation. This is not 
only true for the homogeneous Bianchi type I models but for any homoge- 
neous or inhomogeneous cosmological spacetime which may be cast into the 
"semi-conformal" form (2) with the metric functions depending on z coordi- 
nate as well (technically these are called the spacetimes with two commuting 
orthogonal spacelike Killing vectors and include Bianchi types I, III, V, VI 
and their unidirectional inhomogeneous generalizations). 

Substituting <fi given by the Eq. (22) and the form of the potential into 
the Eq. (18) we are left with a single equation 

GG 2 - GGG + (± - ^ GG 2 + (™ 2 + y) G = 0. (23) 
The substitution 

G = y(G) , (24) 

reduces the degree of the last equation, and after redefining some of the 
constants we finally get 

S+U-f) ^ = 0, (25) 



where 



K=*-\, M = m 2 + ^. (26) 

and y' = dy/dG . 

The Eq. (25) has the following first integral 

Gy' + (K-l)y + — = B, (27) 

y 

where B is an arbitrary constant. Integrating the Eq. (27) one finally obtains 



G = 
where 



2(K - l)y + B - y/A 



B/2(K-1)V~A 



-1/2(^-1) 



{N[{K-l)y 2 + By + M}\ 

(28) 



2(K- l)y + B + y/A_ 

A = B 2 — A(K — l)M , (29) 
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and N is yet another arbitrary constant of integration. The Eq. (28) as it 
stands is valid for A > 0. For A < the integral takes a different form. 

This formally concludes the solution. For, once G(y) is given by the 
Eq. (28) one may in principle invert the expression to get the y(G), then 
using t = J and inverting again to obtain G(t). From G(t) one easily 
reconstructs all the metric functions: p(t) from the Eq. (15), (j)(t) from the 
Eq. (22) and f(t) from the Eq. (12) for example. 

Before turning to analyse the general solution given by the Eq. (28) we 
will present some explicit particular cases of interest in the following Section 
IV to illustrate the above mentioned procedure. 



4 Explicit Exact Solutions 

In this Section we obtain explicitly some exact solutions of interest and briefly 
discuss their behaviour. We start with the simplest ones. 



4.1 The FRW Universes 

The homogeneous and isotropic universes are obtained if one specifies M = 
and B = in the Eq. (27). Note that M = automatically means a = in 
the Eq. (15) which in turn excludes the transversal part of the gravitational 
field p restricting the class of models to LRS ones. If, moreover, the constant 
5 = one always finishes with an isotropic solution. 



If K = (k = 2) one obtains the so called "coasting solution" |10 



A 2 

ds 2 = —e At {-dt 2 + dx 2 + dy 2 + dz 2 ) (30) 

and 

<P = ~At (31) 

This universe expands linearly in synchronous coordinates and has zero de- 
celeration parameter. 

For K ^ one gets from the Eq. (28) 

G = t\ A = ^ (32) 
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Two different classes of solutions appear. If A = 1 one has A = and the 
solution is that described by the Eq. (20) representing the massless minimally 
coupled scalar field FRW universe. If, however, A ^ 1 one gets 

G = t\ ef = ^^t\ 0=-^logt, (33) 

and the metric in synchronous form is given by 

ds 2 = -dT 2 + T 4 / fc2 (dx 2 + dy 2 + dz 2 ) (34) 

The sign of the deceleration parameter for the models described by the 
Eq. (34) depends on the quantity k 2 — 2 so that the model inflates for k 2 < 2 
while decelerates for k 2 > 2. 

4.2 The LRS Models 

Assuming now M = but B ^ and taking K = to obtain 
expression we get from the integral (28) 

G = e* + B, ef = -L e «, <t> = -\ Iog(e* + B). 

The metric in the synchronous form is given by 

ds 2 = -dT 2 + T 2 dz 2 + {T 2 + B) (dx 2 + dy 2 ). (36) 

Note that at early times this LRS model is anisotropic while at T — > oo 
approaches the linearly expanding FRW universe. 
The sign of the deceleration parameter is given by 

- S(e« - f ) , (37) 

and if the constant B < the model never inflates! Note, however, that the 
range of the parameter B might be restricted by the sign of the potential V. 
In the solution (36) the constant B measures the deviations from isotropy 
and as long as B is negative the model does not inflate independently of how 
negligible the deviation from the isotropy is. 

In the case B > the cosmological model decelerates until t < log ^ but 
inflates for later times. Another interesting observation for this model is that 



an analytic 
(35) 
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they are non singular for B > 0. This does not contradict the singularity 
theorems for, if B > one has 

p + 3 p = - { bTt^ <0 ' (38) 

and the energy condition is broken. 

4.3 Non Symmetric Solutions 

We put M^0,5 = 0,Af= — 1/M as well as specify K = again to obtain 
an analytic expression. With these constants the integral (28) gives 

G = sinhwt, u = \fM. (39) 

The transversal and longitudinal degrees of freedom are then given by 

p = — log (tanh — J (40) 

and 

/ = log(smh w*) - — log (tanh — \ + log — (41) 
respectively. And the scalar field evolves according to the following law 

(/> — — log(sinhcjt) + — log ("tanh ^\ . (42) 
2 uj V 2 / 

It is interesting to see from this solution that the anisotropy is contributed 
by both the scalar field (m 7^ 0) and the transversal degree of the gravita- 
tional field (a ^ 0). This can be seen as well from the Eq. (23) where the 
contribution to the anisotropy due to the transversal degree of the gravita- 
tional field and that due to the scalar field enter symmetrically. 

These solutions decelerate for the times t < t c where t c is given by 

1 / ll m 2 + 3 a 2 \ 

t c = - cosh" 1 6 - 4 (43) 
uj \m\j2m 2 + a? ) 

and then, after this time the solutions inflate. 
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5 The Asymptotic Late-Time Behaviour of 
the Generic Solutions 



5.1 The Numerical Analysis 

We now turn to study the asymptotic behaviour of the generic solutions of 
the Eq. (23). From the previous Section we have seen that two different types 
of late-times behaviour occur in exact solutions 

G ~ t N , (44) 

and 

G ~ e Nt . (45) 

These two different asymptotic behaviours correspond to FRW (Kasner if 
N — 1) and to anisotropic models respectively. Technically we have found it 
very difficult to study the integral (28) analytically. In most of the cases one 
can not integrate the Eq. (28) and further integrals in terms of elementary 
functions. We therefore have been forced to use numerical methods. 

To ascertain as to whether the asymptotic behaviour described by the 
Eqs. (44) and (45) occurs in other solutions of the Eq. (23) we have mon- 
itored the values of the following functions while integrating the Eq. (23) 
numerically [fLlp : 

^i(^) £~<2 GG ' 

together with 

n 2 (t) = |. (47) 

It is obvious that if a solution of the Eq. (23) asymptotically follows 
the power law given by the Eq. (44) it is necessary and sufficient that the 
function ri\ defined by the Eq. (46) tends to a constant as the time increases. 
One could use a different function, say n(t) = tG/G, but our numerical 
experiments have shown that the final answer is the same. The quantity in 
(46) was finally retained because many test cases gave the answer somewhat 
faster. Similarly the asymptotic behaviour given by the Eq. (45) is monitored 
by the test function n 2 (or GjG which finally gives the same answer). 



10 



We have integrated numerically the Eq. (23) for different values of the 
constants k, m and a. In each case different initial conditions for G, G and G 
were chosen, though always remaining positive for physical reasons. Negative 
G would imply a change in the signature, negative G would imply a negative 
cosmological constant via the Eq. (12) and finally the positivity of G accounts 
for initial expansion rather than contraction. 

Our numerical results depict a clear scenario for the asymptotic behaviour 
of the Eq. (23). We have found that for every numerical solution the con- 
structed test function n\ tends towards a constant value which depends only 
on the slope of the scalar field potential defined by the constant k. This 
behaviour and the value of the exponent N is absolutely independent of the 
level of the anisotropy introduced either via the scalar field, related to the 
parameter m or purely geometrical anisotropy related with the parameter a. 
For k 2 < 6 the generic solution behaves as 

(48) 

as t goes to infinity which is an asymptotic behaviour of the isotropic FRW 
model. 

For k 2 > 6 the following set of solutions of Eq. (23) 

G = Ct + D, (49) 

which corresponds to an asymptotic behaviour of the vacuum Kasner solu- 
tions or to the solutions with a minimally coupled massless scalar field, is 
an attractor of the Eq. (23) and describes the asymptotic behaviour of its 
generic solution. It is worthwhile to mention that the status of expressions 
(48) and (49) is different: while (48) is only an asymptotic solution, (49) is 
an exact solution for an arbitrary set of constants K, M and B. 

We have also found that the asymptotic behaviour described by the 
Eq. (45) happens only in the case k 2 = 2 and is structurally unstable, for any 
small deviation in the parameter k was changing the asymptotic behaviour 
of the solution to that described by the Eq. (48). 

5.2 The Qualitative Analysis 

The outcome of our numerical calculations can be made plausible by a qual- 
itative analysis of the Eq. (23), which in terms of the new variables 

x = \ogG, y = G, (50) 
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reduces to the following autonomous equation 



y" + (K- l)y' = M- 




y 



(51) 



In the particular case in which M — 0, the resulting linear equation 



If K > 1 (i.e., if k 2 > 6 ), and we assume that G, and thus x, increases 
with time, the asymptotic behaviour of the general solution is y ~ C, which 
corresponds to the set of solutions given by the Eq. (49). 

In the non-linear case, M ^ 0, one could expect that the asymptotic 
behaviour is still the same because if y — > C the right hand side of the 
Eq. (51) goes to 0. Indeed, our numerical experiments show that this is 
the case. This may be seen as well by direct integration of the Eq. (51), as 
depicted in the Fig. 1 for K = 2 and M — 0.1. It is clearly seen that after 
some time the solutions of the Eq. (51) become to a high precision straight 
lines with a slope of -1 (i.e., solutions of the linear equation (52)) and, finally, 
approach the attracting line y' — 0. 

Though unchanging the asymptotic behaviour, the right hand side of the 
Eq. (51) is, of course, important for small y. To stress, it we have drawn as 
dashed lines some solutions of the linear equation (52) corresponding to the 
same initial conditions. 

As for the assumption on G being an increasing function of t, we see in 
the same figure that solutions corresponding to positive initial conditions (as 
required by physical reasons) never reach negative values for y = G, for they 
never cross the line y' — 0. This behaviour was observed in all numerical 
integrations. 

The cases corresponding to k 2 < 6 (—1/2 < K < 1) can be discussed in a 
very similar way by using 



y" + (K - l)y' = 



(52) 



has the following general solution: 



y = C + Ee {l ~ K)x . 



(53) 



x = logG, 
which leads to the following equation 



y = Ge^-V*, 



(54) 



y" + (l-K)y' = M 



2 y' + (l-K)y e2{K - n 

y 2 



(55) 
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Again, the particular case in which M = is trivial and its general solution, 
y = E + Fe^ K ~^ x , will approach y = E, which corresponds to (48). If M ^ 0, 
y = D is no longer a solution of the Eq. (55) for a finite x, but it can still 
represent the asymptotic behaviour of the typical solution, because the right 
hand side of the Eq. (55) vanishes when y — > E and x — > oo. We have 
checked by numerical integration of both the Eqs. (48) and (55) that this is 
indeed the case. 

6 Conclusions 

We have presented here a general exact solution of the Eisntein field equa- 
tions for the anisotropic Bianchi type I universes filled with an exponential- 
potential scalar field. Some of the representative cases and their behaviour 
were explicitely considered in the Section IV. We must stress that by con- 
sidering and studying exact analytic examples one gets a good idea as to 
the behaviour of the general model. Exact solutions, in our view, are inde- 
spensable and must be considered, if possible, before any numerical analysis 
is undertaken to be sure that the results based on numerical "experiments" 
are of any significance. 

Looking at the exact examples we have found three different typical late- 
time behaviours of the models: 

i) The Kasner-like behaviour which is characteristic to the vacuum and mass- 
less scalar field models. 

ii) The FRW-like behaviour, and 

iii) the limiting "coasting" type behaviour. 

Integrating numerically the general solution we have seen that the limiting 
G ~ e Nt behaviour is structurally unstable. This in fact is very interesting, 
for, exactly this type of behaviour as we shall see in paper II, is generic for 
the models with one-dimensional inhomogeneity. So, while in such simple 
models as Bianchi I the instability of these solutions do not cause any reason 
to worry about the isotropization, in more complicated models their stability 
causes problems. 

We have seen that when the constant k 2 , defining the slope of the po- 
tential is less than 6, the generic late-time behaviour is that of an isotropic 
FRW model. In such situations the scalar field acts similarly to a positive 
cosmological constant. 
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For k 2 > 6 the late-time behaviour is that of a Kasner-type universe or, 
which is the same, of the model filled with massless minimally coupled scalar 
field. As long as k 2 > 6 one has no reason to believe that the model will 
isotropize. 

To see whether the models inflate one may look at the sign of the decel- 
eration parameter q. After some algebra and using the first integral given by 
the Eq. (27) the sign of the deceleration parameter q is given by 



1 KG 2 - f km + ^) G + \ (km - B) 2 + \ M. (56) 

2 V 2 / 6 2 



It is easy to see from this expression that as long as the anisotropy parameters 
M and B are switched off (m = a = 0), the inflation of the solutions depends 
only on slope of the potential given by k. Nevertheless, if the anisotropy is 
present the inflation is not only driven by the parameter k but depends as 
well on the rates of anisotropy. Studying the behaviour of the Eq. (50) one 
finds out that solutions generically inflate for k 2 < 2 confirming previous 
results. For k 2 > 2 most of the solutions do not inflate, yet depending on the 
rate of the anisotropy one may find inflating solutions. 
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Figure Caption 



Fig. 1 Phase-space of the Eq. (51) for K = 2 and M — 0.1. The solutions 
correspond to initial conditions in the form y = 0.01, y' = 0.01 + O.ln, with 
n — ... 12. The dashed lines represent the solutions of the linear case (52) 
for the same initial conditions for n — 1 ... 7. 
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